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Peering through the stellar wind of IGR J19140+0951 with 
simultaneous INTEGRAL/ RXTE observations 
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ABSTRACT 

We have used the RXTE and INTEGRAL satellites simultaneously to observe the 
High Mass X-ray binary IGR J19140+0951. The spectra obtained in the 3-80 keV 
range have allowed us to perform a precise spectral analysis of the system along its 
binary orbit. The spectral evolution confirms the supergiant nature of the companion 
star and the neutron star nature of the compact object. Using a simple stellar wind 
model to describe the evolution of the photoelectric absorption, we were able to restrict 
the orbital inclination angle in the range 38°-75°. This analysis leads to a wind mass- 
loss rate from the companion star of ~ 5x 10~ 8 MQ/year, consistent with an OB I 
spectral type. We have detected a soft excess in at least four observations, for the first 
time for this source. Such soft excesses have been reported in several HMXBs in the 
past. We discuss the possible origin of this excess, and suggest, based on its spectral 
properties and occurrences around the superior conjunction, that it may be explained 
as the reprocessing of the X-ray emission originating from the neutron star by the 
surrounding ionised gas. 
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1 INTRODUCTION 

High Mass X-ray Binaries (HMXBs) are binary systems con- 
sisting of a compact object orbiting a massive companion 
star. Prior to the launch of the INTErnational Gamma-ray 
Astrophysics Laboratory {INTEGRAL) in 2002, a large ma- 
jority of the known HMXBs contained a Be companion. In 
Be-type HMXBs, the compact object emits strong X-ray 
flashes when it crosses the equatorial plane of the compan- 
ion star, where a thick disk of matter originating from the 
stellar wind is present. Supergiant O and B-type stars have 
more isotropic stellar winds which absorb the X-ray emis- 
sion of the compact object, rendering them almost unde- 
tectable below a few keV. Thanks to its sensitivity in the soft 
gamma-ray range, however, INTEGRAL has found many 
such systems in the p ast few years (see e.g. iLiu et al1l2006l ; 
iBodaghee et al.ll2007h . In this perspective, the use of X-ray 
spectroscopy at different orbital phases makes it possible to 
probe the stellar wind, providing two-dimensional informa- 
tion on the density and ionization structure of the wind. 
For instance, the soft excess that is present in the soft X-ray 



spectra of many HMXBs, whose origin is still quite mysteri- 
ous, is linked to the physics of the wind close to the compact 
object, especially the region where the fast moving stellar 
wind collides with the slow movi ng and highly ioni zed gas 
surrounding the compact object (jHickox et ajJ|200J ). 

IGR J191 40+0951 was discovered by INTEGRAL in 
March 2003 (|Hannikainen et alj 120031 '). during the first 
observation of the nearby microquasar GRS 191 5 + 1 0^3 
jHannikainen et all l2004h . ICorbet et ail (|2004h and 
IWen et alj (120061 ) analysed archival RXTE / 'All-Sky Moni- 
tor (ASM) data from 1996 onwards and discovered a 13.552 
± 0.006 day periodicity in the X-ray light curve. They 
associated it with the orbital period of the system, identified 
as an X-ray binary. An early analysis of the first RXTE 
and I NTEGRAL data is presented in iRodriguez et al.1 
(2005J, hereafter Paper I). These observations were not 
simultaneous, but the study of the spectral evolution 
of the source allowed them to classify it as a probable 
HMXB hosting a neutron star, which captures matter 
from the stellar wind without Roche lobe overflow. In this 
paper, we confirm this classification for IGR J19140+0951. 
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1 IGR J19140+0951 lies ~ 1.1° from GRS 1915+105. 
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There is no evidence in the RXTE data of an eclipse of 
the compact objet by its companion star, which gives 
an upper limit to the inclination of the system (see Sec. 
4.311. After refin e ment of the X-ray position with Chandra, 
in't Zand et al.l l|2006l) identified the infrared count erpart 
with 2MASS 19140422+0952577. iNespoli et all (|2007f ) later 
estim ated the spectral type o f the companion star to be 
Bl I. lHannikainen et ah! (|2007T ) made further observations, 
narrowing the spectral type to BO. 5 I, confirming the 
HMXB nature of IGR J19140+0951. IGR J19140+0951 
is situated in the direction of the Sagittarius arm, a site 
of active stellar formation. Given that HMXBs are young 
systems and that the Sagittarius arm already hosts several 
HMXBs, it is reasonable to think that IGR J19140+0951's 
distance lies between 2 and 6 kpc. This is further supported 
by the upper l imit of 5 kpc obtained fro m the K-magnitudes 
of the source (|Hannikainen et al.l l2007). 

In this paper, we present several RXTE and INTE- 
GRAL observations of IGR J19140+0951, most of them 
taken simultaneously with the two satellites. The sequence 
of observations and data reduction process is described in 
Section [2] In Section [3] we present the spectral analysis and 
its main results. In Section |4j we discuss these results and 
use a simple stellar wind model to calculate some character- 
istics of the source. 



2 OBSERVATIONS AND DATA REDUCTION 

We have been monitoring GRS 1915+105 (and therefore 
IGR J19140+0951) since March 2003 with INTEGRAL, and 
made simultaneous RXTE observations as often as possible. 
Our data cover the period March 2004-September 2007 with 
32 observations, the log of which is reported in Table [T] We 
also used one public RXTE ob servation performed in March 
2003, and used the results from lin't Zand et alj l|2006l ) based 
on one Chandra observation taken in May 2004. 



2.1 INTEGRAL/IBIS data 

To study the properties of the hard (> 20 keV) X-ray emis- 
sion from the source, we use d the data from the IBIS tele- 
scope on board INTEGRAL l|Ubertini et al.ll2003T ). This in- 
strument uses a coded mask, which allows imaging over a 
large Field of View (FOV), ~ 30° x30° up to zero response. 
The Totally Coded FOV (TCFOV), where the sensibility is 
uniform, is 9 °x9°. The INTEGR AL Soft Gamma-Ray Im- 
ager (1SGRI, iLebrun et al.l [20031 ) is the upper layer of the 
IBIS detection unit, covering the range between 13 keV and 
a few hundred keV. Its angular resolution is 12'. 

The INTEGRAL campaign was aimed at studying the 
microquasar GRS 1915+105, but also allows observations 
of every other source that lies in the IBIS TCFOV, includ- 
ing IGR J19140+0951. Since October 2004, our RXTE ob- 
servations were made simultaneously with the INTEGRAL 
observations. Following the method used in Paper I, the 
data were reduced using the standard Off-line Scientific 
Analysis (0SA) v. 7.0 software package provided by the IN- 
TEGRAL Science Data Centre ( |http://isdc. unige.ch). First 
we ran the software up to the production of images and 
mosaics in the 20-40 and 40-80 keV range, using only 



the science windows simultaneous with the RXTE observa- 
tions. The software was left free to find the most significant 
sources in the images, which means that faint sources like 
IGR J19140+0951 were not always detected. Then we used 
a catalogue containing the 7 most luminous sources of the 
field for spectral extraction, in order to "force" the program 
to extract the spectrum of our source from the images. Given 
the faintness of the source, we rebinned the redistribution 
matrix file provided with OSA to obtain 6 energy bins in 
the range 18-110 keV. 

2.2 B1T£/PCA data 

For our analysis of the RXTE data, we used the top layers 
of the Proportional Counter Units (PCUs) 2 and 3 of the 
Proportional Counter Array (PC A), which were turned on 
in every observation. The data were reduced using LHEAS0FT 
package v6.4. We restricted the study to the time when the 
elevation angle was above 10° and the offset pointing was 
less than 0?02, as recommended by the GOF (Guest Ob- 
server Facility) to avoid Earth contamination and errors due 
to slew motion. In addition, and since the source is quite 
weak, we further rejected times of high electron background 
in the PCA (i.e. times when the electron ratio in one of the 
PCUs is greater than 0.1). The response and background 
were generated using PCARSP version 10.1, and PCABACKEST 
version 3.1 respectively, with the latest calibration files pro- 
vided by the RXTE GOF. These files include the corrections 
applied to the PCA Background since 2007 September 18. 
We verified for IGR J19140+0951 that the impact of these 
corrections on the background is typically 5% below 5 keV, 
and 2% above, which is what is expected for a faint source. 
The so urce is faint and th e errors dominated by statistical 
effects (|jahoda et al.|[2006f ). so we did not include any sys- 
tematic error. We however checked that the results remained 
the same when adding a "standard" systematic error of 0.6% 
to the data, with only a slight and non significant decrease 
in the reduced \ 2 . in addition, the default spectral binning 
of the PCA detector is too fine for this dim source, so we 
rebinned all our RXTE /PCA spectra above 8 keV to make 
one bin out of 4 native bins between 8.1 and 16.5 keV, and 
one bin out of 10 bins between 16.5 and 26.6 keV. This im- 
proves slightly the precision on the parameters. 

We also checked the data from the High Energy Tim- 
ing Experiment (HEXTE), the second instrument on-board 
RXTE. However, as already pointed out in Paper I, no 
HEXTE data can sensibly be used in our analysis because 
of the faintness of the source. 

The PCABACKEST software computes the to- 
tal instrument and cosmic X-ray background. Since 
IGR J19140+0951 is a faint source, and lies in the Galactic 
Ridge, we also correc ted the spectra for t he Ga lactic X-ray 
background (GXB). IValinia fc Marshalll (|l998l ) quantified 
the diffuse GXB in several regions of the Galactic ridge, 
including one with Galactic Latitude in the range -1.5-1.5° 
and Galactic Longitude in the range -40-40°, which corre- 
sponds to our source. They used a three component model 
in the 3-35 keV range: a Raymond-Smith model and a 
power law, attenuated by interstellar absorption with a col- 
umn density of Nh 1s ~ 1.5 10 22 cm -2 . Note that the model 
absorption measures the Galactic column density through 
the whole Galaxy in a given direction, so it overestimates 
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Table 1. List of the observations. 



Observation ID 


INTEGRAL 


Date_ 


RATE Exposure 


Orbital Phase_ 




revolution 


(MJD) 


(s) 


(cycles) 


P80404-01-01-00 


0049 


52,708.847 


2864 


0.815 


P90112-01-01-00 




53,087.593 


7760 


0.763 


P90112-01-02-00 


0246 


53,296.931 


8752 


0.210 


P90112-01-03-00 


0295 


53,443.185 


1936 


0.002 


P90112-02-01-00 


0305 


53,473.889 


5808 


0.268 


P90112-02-01-01 


0305 


53,473.202 


3232 


0.217 


P90112-02-02-00 


0361 


53,641.355 


2016 


0.625 


P90112-02-03-00 




53,676.250 


3216 


0.199 


P90112-03-01-00 


0373 


53,677.239 


9152 


0.274 


P90112-03-02-00 


0431 


53,852.034 


2352 


0.170 


P91083-01-01-00 


0478 


53,990.140 


1776 


0.361 


P91083-01-01-01 


0478 


53,990.207 


2176 


0.366 


P91083-01-02-00 


0480 


53,996.118 


2416 


0.803 


P91083-01-03-00 


0481 


54,001.625 


2272 


0.209 


P91083-01-04-00 


0486 


54,016.558 


1488 


0.311 


P91083-01-05-00 


0489 


54,025.374 


2176 


0.961 


P91083-01-06-00 


0493 


54,037.379 


1408 


0.847 


P91083-01-06-01 


0495 


54,040.956 


1520 


0.111 


P91083-01-07-00 




54,049.996 


4816 


0.778 


P91083-02-01-00 


0501 


54,059.039 


3200 


0.445 


P91083-02-02-00 


0537 


54,166.928 


1680 


0.406 


P91083-02-03-00 


0539 


54,173.148 


1520 


0.865 


P9 1083-02-04-00 


0542 


54,181.948 


4336 


0.514 


P91083-02-05-00 


0544 


54,187.806 


6064 


0.947 


P91083-02-06-00 


0546 


54,195.604 


1904 


0.523 


P91083-02-06-01 


0546 


54,195.537 


1504 


0.518 


P91083-03-01-00 


0549 


54,202.729 


6832 


0.048 


P9 1083-03-02-00 


0551 


54,208.841 


3216 


0.499 


P91083-03-03-00 


0553 


54,214.537 


1920 


0.920 


P91083-03-04-00 


0557 


54,226.594 


3040 


0.809 


P9 1083-03-05-00 


0559 


54,232.526 


2160 


0.247 


P91083-03-05-01 


0559 


54,232.594 


2528 


0.252 


P91083-04-01-01 




54,373.323 


2720 


0.636 


Chandra 




53,136.732 




0.387 



a Midpoint of observations 

k The orbital phase s are calculated using an orbital period of 13.552 days, and a phase 0.5 in MJD 51593.4 
llCorbet et al.ll2004l) 



the absorption experienced by IGR J19140+0951 which is 
most probably situated between 2 and 6 kpc. As we are 
interested in the relative variations of the absorption, and 
since the Galactic absorption in the direction of the source 
is constant, the precise value of the latter is not needed as 
it just adds a constant to the absorption intrinsic to the 
source. 



Moreover, the smallest measurement of the absorption 
we have near phase with a meaningful error bar is 3.9 ± 
2.3 10 22 cm -2 , once the Galactic contribution has been re- 
moved. This means that, even if we assume the minimum 
intrinsic absorption is ~ 1.6 10 22 cm -2 , this is still bigger 
than the Galactic absorption in the direction of the source. 
So we can reasonably consider the absorption we measure 
as highly dominated by intrinsic absorption. In the follow- 
ing, we assume that what we measure is the local absorp- 
tion, which is, hereafter, referred to as Nu- The RXTE /PCA 
spectra were fitted in XSPEC version 11.3.2, between 3 and 



25 keV. A joint RXTE-INTEGRAL spectrum is shown in 
Fig.CD right. 

3 RESULTS 

3.1 Spectral analysis 

We used the ephemeris of ICorbet et all (|2004T ) to fold the 
spectra, starting at phase when the flux is at a minimum. 
The analysis of the absorption underwent by the X-ray flux 
along the orbit shows that phase also corresponds to a 
minimum of the absorption (see section [4.3[) . Thus, phase 
corresponds to when the compact object is located be- 
tween the Earth and the companion star (hereafter inferior 
conjunction): in this case, the X-ray flux coming from the 
compact object travels a shorter distance in the compan- 
ion stellar wind. Phase 0.5 is therefore when the compact 
object is behind the companion star (superior conjunction). 
The fact that a minimum flux corresponds to a minimum 
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Figure 1. Left: Residuals in terms of <r between the models used to fit the INTEGRAL+RXTE spectra. In each case the model 
mentioned in the panel is convolved with interstellar absorption. Po stands for power-law, high stands for high-energy cut-off, gaus stands 
for Gaussian and bb stands for black body. Right: Joint INTEGRAL+RXTE spectrum from ObsID P91083-03-02-00, corresponding 
to the residuals on the left. RXTE/PCA data are from 3-25 keV and INTEGRAL /ISGKl data from 20-100 keV. The best model is 
superposed on the spectrum. 



absorption may be puzzling, but actually these two effects 
act on two different energy ranges: indeed, the mo dulation of 
the A SM flux is seen only in the 5-12 keV range ()Wen et al.l 
2006), where the photoelectric absorption has no influence. 
Qualitatively, when the neutron star is behind the compan- 
ion star, we see directly the shock between the stellar wind 
and the ionized gaz surrounding the compact object, where 
the high-energy emission is produced. At inferior conjunc- 
tion, however, the high energy emitting region is partially 
occulted by the neutron star: thus, this moment appears as 
a minimum flux in the ASM. 

The phases are calculated at the mid-point of the ob- 
servations. One observation typically lasts 0.005 phase. The 
0.003 d ay uncertainty in the period determination of 13.552 
days bv lWen et al.l l|2006l ) leads to an uncertainty of 0.02 on 
the phase determination over the 4 years of observations. 

We tested several models in analysing the spectra, start- 
ing with phenomenological models. The data were well fitted 
using an absorbed power law combined with a high energy 
cut-off. An iron fluorescent line at 6.4 keV is present in the 
observations where the source is detected at a high enough 
significance (see Fig. [1] right). Note that the high energy 
cut-off is only detected in the observations where the source 
flux is sufficiently high, thus giving enough precision at the 
high energy end of the spectra. Fig. (TJ left, shows the resid- 
uals from the fit procedure when we add successively each 
component of the model. Note that for this observation, we 
need an additional component to accurately describe the 
spectrum below 6 keV. The spectrum shows a "soft excess" 
feature in this range, which is discussed below. 



3.2 Absorption evolution 

The most important feature that we found is a strong evolu- 
tion of Ah along the orbit (Fig. [2}. In order to check whether 
the evolution of Ah was genuine or not, we produced power 
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Figure 2. Evolution of the absorbing column density along the 
orbital phase of the system. Error bars are shown at the 90% con- 
fid ence level. "Chandra d ata" refers to the value of jVjj measured 
bv lin't Zand et al.l ll2006h . 



law photon index (r) vs. Ah contour plots for all obser- 
vations. Eight typical observations are reported in Fig. [3] 
These correspond to spectra whose parameters are well con- 
strained. The outer contour represents the 99% confidence 
level. The parameters T and A^h are correlated: any vari- 
ations of Ah will have a strong impact on the measure of 
F and vice-versa. The elliptical elongation of each individ- 
ual contour is the result of this link in the fitting process. 
The contours are clearly distinct from one observation to the 
next, which clearly indicates that the evolution of A^h and 
F are genuine. If this were not the case, the contours would 
overlap. The detailed study of the evolution of Ah along the 
orbital phase is presented in Sec. 14.31 

The absorption measured by Chandra is more precisely 
determined because of the better spectral resolution and the 
lower energy boundary of this instrument in the soft X-ray 
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Table 2. 



Table 2. Best fit parameters. The errors and upper limits are given at the 90% confidence level. * observations where a soft excess is 
detected with a 90% confidence level; **with a 99.9% confidence level. The reduced x 2 are calculated for 16 degrees of freedom (d.o.f.) for 
observations with both RXTE and INTEGRAL, and 11 d.o.f. otherwise. The unabsorbed flux is calculated in the range 2-20 keV. c marks 
the Chandra data and R marks the 4 observations for which INTEGRAL data are missing. 



Phase 




r 


E c „i 


Unabs. flux 


xi 


N u 




Tp 


y 




xi 




(10 22 cm- 2 ) 




(keV) (10" 


11 erg cm -2 s -1 ) 




(10 22 cm- 2 ) 


(keV) 


(keV) 


^ 4feT e 


r 2 ) 
























Phenomcnological model 








Comptonization model 










0.002 


< 9.8 


2.2 ± 0.37 




2.6 ± 0.6 


1.79 
















0.048 


< 2 


1.68 ± 0.18 


7.5 ± 2.1 


6.2 ± 2.2 


1.19 
















0.111 


< 5.9 


2.45 ± 0.29 




4.2 ± 1 


1.19 
















0.170 


3.9 ± 2.3 


1.84 ± 0.16 




9.5 ± 1 


1.42 
















0.199 H 


6.9 ± 2.1 


2.47 ± 0.17 




6.5 ± 0.5 


1.20 
















0.209 


< 13.2 


2.35 ± 0.75 




0.8 ± 0.1 


1.00 
















0.210 


< 1.6 


1.25 ± 0.3 


6.4 ± 0.6 


3 ± 2 


1.39 


< 4.5 


3.7 ± 0.7 


6.4 ± 1.5 


1.1 ± 


0.6 


1 


23 


0.217 


< 4.6 


1.05 ± 0.36 


6.2 ± 0.7 


9.2 ± 0.8 


1.01 


6.1 ± 2.0 


3.6 ± 0.8 


5.9 ± 1.2 


1.0 ± 


0.5 


1 


63 


0.247 


3.8 ± 1.6 


1.35 ± 0.15 


6.3 ± 0.5 


19 ± 5 


1.80 
















0.252 


< 4.8 


1.60 ± 0.35 


< 18 


5.2 ± 2 


1.61 


< 4.0 


5.5 ± 1.6 


4.6 ± 1.3 


0.8 ± 


0.5 


1 


56 


0.268 


12.4 ± 1.7 


1.15 ± 0.1 


9.2 ± 0.6 


67 ± 3 


1.08 


34.3 ± 1.8 ** 


13.3 ± 2.6 


2.0 ± 0.4 


0.4 ± 


0.2 


1 


12 


0.274 


6.3 ± 1.4 


2.3 ± 0.1 




12 ± 2 


1.29 
















0.311 


11.2 ± 3.1 


2.15 ± 0.19 




7.6 ± 1 


1.13 
















0.366 


18.7 ± 2.0 


1.57 ± 0.15 


6.1 ± 0.8 


27 ± 2 


1.62 


27.1 ± 5.0 ** 


13. ± 4 


2.6 ± 0.5 


0.6 ± 


0.3 


1 


00 


0.387 c 


10. ± 0.3 


1.1 ± 0.8 




35 ± 5 


1.01 
















0.407 


18.2 ± 3.0 


2.03 ± 0.14 




12 ± 1 


1.27 


17. ± 3.0 


11.5 ± 7.0 


< 4.66 


0.4 ± 


0.8 


1 


30 


0.445 


17.6 ± 8.0 * 


1.0 ± 0.6 


< 10.4 


32 ± 1.5 


1.18 


32. ± 4.5 * 


10.7 ± 4.1 


3.6 ± 0.8 


0.7 ± 


0.5 


1 


04 


0.499 


21. ± 10.** 


0.9 ± 0.15 


12.2 ± 0.8 


13 ± 2.5 


1.27 


< 26.6 


5.1 ± 0.5 


9.5 ± 4.5 


3.0 ± 


3.0 


1 


41 


0.514 


24.5 ± 3.9 


1.45 ± 0.27 


9.3 ± 1. 


14 ± 3 


1.15 


28.3 ± 2.5 


4.9 ± 0.3 


5.5 ± 1.0 


0.9 ± 


0.4 


1 


74 


0.518 


28.4 ± 2.1 


1.93 ± 0.12 


10.2 ± 2 


48 ± 5 


1.60 


25. ± 2.2 


5.5 ± 0.7 


4.8 ± 0.8 


0.9 ± 


0.3 


1 


66 


0.523 


22.4 ± 4.3 


2.02 ± 0.2 


< 30 


15 ± 4 


1.34 


21.6 ± 3.8 


15 ± 10 


2.2 ± 1.2 


0.5 ± 


0.8 


1 


30 


0.625 


10.2 ± 1.8 * 


1.6 ± 0.2 


6.6 ± 0.9 


27 ± 5 


1.16 


13.4 ± 1.0 ** 


13 ± 5.0 


2.5 ± 0.7 


0.5 ± 


0.4 


1 


87 


0.636 fl 


13.1 ± 1.6 


1.14 ± 0.13 


6.3 ± 0.5 


47 ± 5 


1.38 


14.2 ± 1.5 


4.9 ± 0.4 


6.0 ± 0.6 


1.4 ± 


0.3 


1 


04 


0.763 fl 


7.9 ± 1.7 


1.16 ± 0.2 


6.2 ± 0.5 


16 ± 3 


1.67 


9.1 ± 1.5 


4.3 ± 0.3 


6.0 ± 0.7 


1.2 ± 


0.3 


1 


05 


0.778 H 


< 3.2 


1.35 ± 0.55 


< 9.2 


6 ± 2 


1.75 


3.2 ± 2.2 


4.9 ± 0.8 


5.0 ± 1.2 


0.9 ± 


0.5 


1 


95 


0.803 


6.0 ± 5.0 


1.25 ± 0.45 


6.2 ± 1.5 


5.6 ± 1 


1.67 
















0.809 


8.7 ± 2.5 


2.32 ± 0.18 




5.8 ± 1 


1.78 
















0.815 


13.2 ± 2.5 


1.95 ± 0.14 




13.2 ± 2 


1.68 
















0.847 


10.4 ± 1.6 


2.28 ± 0.11 




18 ± 2 


1.12 
















0.866 


7.6 ± 2.9 


1.25 ± 0.37 


6.4 ± 0.8 


18 ± 2 


1.48 


8.7 ± 2.4 


4.3 ± 1.0 


5.6 ± 1.2 


1.1 ± 


0.5 


1 


49 


0.920 


< 1.7 


1.52 ± 0.18 


6.2 ± 0.7 


9.6 ± 2 


1.27 


2.3 ± 1.7 


4.0 ± 0.9 


5.6 ± 0.9 


0.8 ± 


0.4 


1 


79 


0.947 


< 20 


2.3 ± 0.5 




0.3 ± 0.2 


1.75 
















0.961 


< 10 


2.4 ± 0.7 




1.1 ± 0.5 


1.23 
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Figure 3. Contour plots of x 2 m the Nn-T plane from fits of eight spectra. The contours correspond to the 68, 90 and 99% confidence 
limit, the central cross being the most probable value. The numbers next to the contours are the orbital phase of each observation. Left: 
Five observations with different 7Vh values and similar F Right: Four observations with similar TVh values and different F 



range. It is ~ 40% lower than the general tendency observed 
by RXTE around the same phase. The discrepancy may 
have two origins. Firstly, there may be a difference of cali- 
bration between the two instruments. Secondly, the Chandra 
high energy range is limited to 10 keV, so it gives bad con- 
straints i_on_the_index_i of ^_the power law, as can be seen in 
Fig. Ill (|in't Zand et all (|2006l ) obtained 7V H =T0 ± 0.3 10 23 
cm -2 and r=l.l ± 0.8). Since increasing V would increase 
JVh, confidence levels in the A^h-F plane could increase the 
value of the error on Ah- Given these uncertainties the value 
of A^h obtained with the Chandra observation was excluded 
from our wind model fits. 

The absorption measured using the Comptonization 
model differs from the phenomenological model for some ob- 
servations, especially in the observations where a soft excess 
is detected, but its overall evolution remains the same. It 
confirms the reality of this evolution, and the discrepancy is 
inherent to the differences between the models. Because the 
phenomenological model fits all the spectra, and in a view to 
obtain a consistent analysis, we used only the values of A^h 
obtained from the phenomenological model in our studies. 

3.3 Thermal Comptonization 

Since a cut-off power law is usually interpreted as a sig- 
nature of thermal Comptonization, we replaced this simpler 
model with a mo re physical Comptonization model (C0MPTT, 
lTitarchukiri994l 'l. In ICR J19140+0951, no stable accretion 
disc nor thermal component has been detected, so we fixed 
the disk temperature parameter kTo at 0.1 keV. The xl ob- 
tained when fitting with this model (modified by photoelec- 
tric absorption) were comparable to those of the phenomeno- 
logical models. Since the parameters are constrained by the 
high energy end of the spectra, only observations with a rel- 
atively high flux, hence a relatively good signal/noise ratio 
at high energies, are suitable for the Comptonization model. 

3.4 Soft Excess detection 

In some observations, the xt was f ar higher than for the 
other observations, the discrepancy between the model and 
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Figure 4. Evolution of the photon index along the orbital 
phase of the system. Error bars are shown at the 90% confi- 
de nce level. "Chandra d ata" refers to the value of T measured 
bv lin't Zand et al.l ||2006T) . 

the data being particularly important at low energies. In 
those cases, the residuals between the model and the spec- 
trum show that some soft excess is present (see Fig.[T] left). 
This featur e has already been observed in many X-ray bi- 
naries (see iHickox et al.l |2004| . for a review) . We tried to 
add systematically a black body component at low energies 
to all our observations, obtaining a better x 2 for several of 
them. In four observations, an F-statistic test indicated val- 
ues greater than 99.9% for a true improvement of the fit with 
the black body, giving a very good confidence in this detec- 
tion. In one more observation, the test gave a value greater 
than 90%. The spectral characteristics of these observations 
are given in Table [3] 

Note that other models, such as a Raymond-Smith 
model, could also account for the soft excess, but we lack 
the resolution at low energies to discriminate between the 
different possibilities, so we kept the simplest one, i.e. the 
black body model. Besides, since the black body and the 
absorption models influence the spectra mostly below 5 keV 
where the PCA coverage and resolution are moderate (only 
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Table 3. Spectral characteristics of the five soft excess observa- 
tions. 



Phase 


kT 




r 






0.268 


0.27 ± 0.02 


12.4 


± 1.7 


1.15 ± 0.1 


1.08 


0.4 


0.366 


0.27 ± 0.10 


18.7 


± 2.0 


1.57 ± 0.2 


1.62 


0.6 


0.445 


0.22 ± 0.04 


17.6 


± 8.0 


1.00 ± 0.6 


1.18 


0.7 


0.499 


0.46 ± 0.25 


21. 


± 10. 


0.90 ± 0.2 


1.27 


3.0 


0.625 


0.27 ± 0.08 


10.2 


± 1.8 


1.60 ± 0.2 


1.16 


0.5 



6 x 10 22 cm -2 , while it is clear from the present study that 
./Vh is tightly linked to the orbital phase. Secondly, the pro- 
cedure of summing spectra taken at similar fluxes in Paper 
1 may have allowed us to pick up peculiar phases of accre- 
tion, as e.g. flares due to enhanced accretion. The analysis 
reported here focuses more on a global and "normal" be- 
haviour of IGR J19140+0951. Finally, since the analysis of 
the INTEGRAL data in Paper I was based on the detec- 
tion by ISGRI, it is clearly biased towards states that are 
bright at hard X-rays. Our orbital-phase dependent analy- 
sis showed us that there was not necessarily a connection 
between the level of hard X-rays and the position on the 
orbital phase. Note also that although IGR J19140+0951 
and GRS 1915+105 are 1.1° apart, some confusion is still 
possible in the JEM-X data. 

Two of the RXTE observations presented here were al- 
ready analyzed in Paper I. They respectively correspond to 
phase 0.763 (Obs. 3 in Paper I) and 0.815 (Obs. 2 in Pa- 
per I). One can immediately see that the results obtained in 
both analyses clearly differ, in particular in the value of the 
absorption (and consequently the unabsorbed fluxes). There 
are two reasons that are the very probable origin of the dis- 
crepancies. First, the analysis presented in Paper I was made 
with old versions of the software and history file that could 
cause errors on the background estimation of up to 10%. 
This problem has been solved since then. The second origin, 
which is probably the main reason for the differences, is that 
in Paper I the Galactic background was not taken into ac- 
count in correcting the data. This has a strong effect at low 
energies, and hence influences highly on the determination 
of the absorption. 

Qualitatively, however, the spectral behaviour of 
IGR J19140+0951 is confirmed. We confirm that an iron 
Kq fluorescence line is present in most of the spectra. This 
line is intrinsic to the source, since if it were due to the 
Galactic background, we would expect its flux to be roughly 
constant, which is not the case here. We measure upper lim- 
its on its width in the range 300-1000 eV, depending on the 
line flux. This could be indicative of a narrow line, rather 
than a broad line, but actually we are mainly limited by the 
instrumental spectral resolution. 

4.2 Spectral evolution 

4-2.1 Comptonization model 

Our first spectral analysis showed that a simple cut-off power 
law model could accurately fit all the spectra. Therefore we 
applied a thermal Comptonization model, as this is the usual 
physical model which explains this feature. We obtained the 
parameters kT e and r reported in Table [5] which are con- 
sistent with similar observations of other HMXBs and ab- 
sorb ed_^mirces_det^tedby INTEGRAL (e.g. IGR J16320- 
475. lRodriguez et al.ll2006h . The Comptonization parameter 
y — — — §-t (also known as the Kompaneets parameter) de- 
termines the efficiency of the Co mptonization pro cess, and 
thus the shape of the spectrum l|Titarchukl Il994r i. In our 
case, this parameter is ~ 0.5-1, indicating a moderately 
efficient process. It corresponds to a rather low accretion 
rate, which is consistent with the supposed nature of the 
system: wind-fed accretion is less efficient than Roche lobe 
overflow. Our estimation of this mass-loss rate (Sec. 14.3 . 2 (I 



a Phenomenological model. 

" Calculated for 15 d.o.f. for the phenomenological model. The 
Comptonization model gives x 2 values with differences of ± 0.15. 
c Kompaneets parameter obtained from the Comptonization 
model. 

6 spectral channels cover the range 3-6 keV), it is difficult to 
separate their respective contributions, and hence we obtain 
large uncertainties for Ah ■ Similarly, the black body temper- 
ature that we found in the range 0.3-0.6 keV must be treated 
with caution. Other X-ray binaries have soft excess temper- 
atures in the range 0.1-0.2 keV; a black body with this tem- 
perature has weak influence above 3 keV where RXTE is 
sensitive, so the values we obtain should be regarded more 
as an order of magnitude than absolute values. 

Adding a black body at low energy increases the un- 
certainties as it adds a degeneracy in the fitting process: 
increasing the black body flux compensates for a decrease 
in the absorption level. Therefore we verified, by drawing 
confidence levels in the kT-A*H plane, that the black body 
temperature is clearly determined. Given the high uncer- 
tainties on this latter parameter, however, the contours of 
several observations in which the soft excess is detected over- 
lap. We, therefore, cannot conclude on any evolution and/or 
relation of the temperature with orbital phase. 



4 DISCUSSION 

4.1 Comparison with previous work 

A comprehensive analysis of the very first X-ray observa- 
tions of IGR J19140+0951 is reported in Paper I. The goal 
of this study was to try to understand the nature of this, 
then newly discovered source. At that time only a set of 
three RXTE observations and a limited number of simul- 
taneous JEM-X and ISGRI data were available. Instead of 
performing an orbital phase resolved study, and in order to 
improve the signal to noise ratio of the INTEGRAL data, 
we performed a time-resolved spectroscopy based on the 20- 
40 keV luminosity level of the source. As a result it is very 
likely that we mixed the spectra extracted at different or- 
bital positions. This is particularly obvious if one tries to 
compare the results obtained in paper I and the ones we re- 
port here. Although qualitatively they match quite well, i.e. 
the spectral models are similar, and the sets of parameters 
are compatible, none of the "states" identified through the 
INTEGRAL data in Paper I were observed here. 

First of all it has to be noted that during the fitting 
of the INTEGRAL spec tra in Paper I, the value of Ah was 
frozen to that found in I Swank and Markwardt (2003]), i.e. 
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in the range 4-7xl0~ 8 MQ/year corroborates the B-type 
nature of the companion. Note that it has no significant 
evolution along the orbital phase, showing that the emission 
process is remarkably stable along the orbit. This suggests 
a rather stable accretion rate, and therefore would tend to 
favour a circular orbit for the system. 

4-2.2 Phenomenological model 

Several observations show a photon index V ~ 1, for instance 
1.15 ± 0.1 at phases 0.268 and 0.636, or 1.05 ± 0.36 at phase 
0.217. This is consistent with a neutron star nature of the 
compact object, as these hard spectra are typical of satu- 
rated Compt onization found in bin ary systems containing a 
neutron star (|Di Salvo et al.ll2004T ). 

Having identified the emission process as thermal 
Comptonization, we returned to the phenomenological 
model, which allows all the spectra to be fitted and permits 
a true study of the absorption evolution to be done. Before 
studying the evolution of Nh, we can draw some conclusions 
from the other parameters. 

The unabsorbed flux undergoes strong variations that 
are not correlated with the orbital phase. The variations 
probably arise from the details of the accretion process. The 
stellar wind is not perfectly homogeneous and the accretion 
onto the neutron star, following the magnetic field line to 
the magnetic p oles, is a highly un stable and complicated 
process (see e.g. iTakata et al.ll200q . and references therein, 
for a review of the emission processes). Since our observa- 
tions are spread over several years, it is not surprising that 
the source flux is highly variable. However, if the orbit of 
the neutron star were eccentric, when the compact object 
approaches the companion star where the stellar wind den- 
sity is higher, we would see an increase in the source flux. 
Such evolution is marginally visible, the unabsorbed flux be- 
ing smaller between phases 0.9 and 0.2. However as this is 
only marginal, we interpret it as another point in favour of 
a circular orbit for the system. 

The photon index of the power-law component under- 
goes strong variations as well (Fig. [3} . These variations are 
not correlated with the flux evolution or the orbital phase, 
which means that this is probably just another effect of the 
high variability of the accretion process. 



4-2.3 Origin of the Soft Excess 

During our study, some spectra exhibited an excess in the 
soft X-ray part of the spectra, which we modele d by adding a 
black body component to the model. Following iHickox et al.l 
|2004) ■ we can look for possible mechanisms explaining this 
feature. Using the parameters listed in Table [4] the orbital 

" G M P 2 1 ^' 3 

separation between the two objects a (a = — 4 7V 'i° rb ) i s 
in the range 51-74 Rq. If we assume a mass for the neutron 
star of ~ 2 Mq, then the Lagrange point Li is located at 
a distance from the companion star in the range 34-56 Rq. 
Thus, it is very unlikely that the companion star overflows 
its Roche lobe, which is confirmed by the fact that we do 
not observe any accretion disk. This rules out explanations 
for the soft excess feature involving an accretion disk or 
a gas s tream flowing from the c ompanion to the compact 
object |Fahlman fc Walker|[l980l ). We can also rule out soft 



emission from the accretion column, as it would imply soft 
pulses in the spec tra. 

According to IHickox et all (|2004r i, this leaves two pos- 
sible processes which involve the immediate surroundings of 
the compact object, characterized by the influence of the 
compact object X-ray emission. The stellar wind of mas- 
sive stars is mostly accelerated by bound-bound transitions 
of atoms. The neutron star, due to its high energy emis- 
sion, ionizes the surrounding material, so the already ion- 
ized gas around the neutr on star is no longer acceler ated by 
the stellar radiation field |Fransson fc Fabianlll980h . When 
the compact object moves along its orbit, the hot gas will 
gradually be overtaken by the stellar wind. This will lead to 
the formation of a "tail" trailing the neutron star. 

The stellar wind collides with this tail, transforming 
a fraction of its kinetic energy into X-ray emission. This 
results in a shock located b etween the compact object and 
the secondary. According to IHickox et alj (|2004j ) . this could 
explain the soft excess of faint sources (Lx < 10 36 erg s _1 ), 
which is the case of ICR J19140+0951 {L x ~ 3.10 35 erg s~ x 
for a distance of ~ 5 kpc and a characteristic X-ray flux of 
~ 10 erg cm -2 s _1 ). Finally, the last process involves the 
"tail" itself: the diffuse cloud around the neutron star may 
reprocess the hard X-rays coming from the compact object. 

However, to differentiate between emission produced in 
collisionally or photoionized plasma, high resolution grat- 
ing spectra are required, which is unfortunately beyond the 
possibilities of RXTE. Still, the fact that the soft excess 
was only seen around the superior conjunction may be of 
some importance. Until now, only a handful of soft excess 
detections have been reported and among them two d etec- 
tions, reported for 4 U 1700-37 bvlHaberl et al.l (|l989T ) and 
for Centaurus X-3 bv lBurderi et all (|2000h 7 both sources for 
which the orbital position is known, fall during phase 0.5- 
0.75. However, this is a too small sample to draw any definite 
conclusion. 



4.3 Wind model 

4-3.1 Description of the model 

Since the orbit of the X-ray source is probably almost circu- 
lar, one possible explanation for the observed change in JVh 
relies on an inclined orbit for the system. Qualitatively, when 
the neutron star is behind its companion its light has to 
travel a longer distance in the stellar wind, so the absorption 
in the soft X-rays increases. Therefore, at the superior con- 
junction we should observe a maximum in the absorption, 
which is indeed the case (Fig. [3]). Following this hypothesis, 
we used a simple wind model in order to constrain several 
paramet ers for the s y stem, following the method already 
used bv lLevine et all (|2004h . A B-type star emits a strong 
stellar wind, us ually taken to be st ationary and spherically 
symmetric (e.g. ICastor et alj|l975h . Its velocity is given by 
a /3-law: 



V(r) =V + (Voo - V ) (l - 



(1) 



where r is the distance from the companion star, 7?* 
its radius, Vq the wind velocity at the photosphere, and 
Voa the terminal wind velocity. Observations show that Voo 
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Figure 5. Model absorption (continuous line) and experimental absorption (A symbols), for the stellar model with parameters M* = 
20Mq , R* = 21Rq and /3 = 0.5. Left :, result obtained with data and model normalised (jq = 1.440). Right : Same as left, but 
without normalisation (xi = 1.432). 




\\\ 






The equation of mass conservation then gives 



= 47ir V(r)p(r) 



(3) 



where M* is the mass-loss rate of the star, and p(r) the 
density. The number density n(r) of hydrogen atoms is in 
turn given by p(r) = fim,HTl(r), where p, is the mean atomic 
weight of the particles expressed in units of the hydrogen 
atom mass mn {p = 0.602 for solar composition material). 
The instantaneous equivalent absorbing column density of 
hydrogen Nh between the neutron star and the observer is 
given by: 



N H = N„ IS 



+ 



N H IS M + 



n[r(s)]ds 



ds 



- 2 (1- 



2*.)t> 

r ' 



(4) 



where Nh ism is the contribution from the interstellar 
medium, and s is the distance along the line joining the 
neutron star and the observer, s is given by 



Figure 6. Most probable orbital inclination, as a function of the 
companion radius, for various M* and p. a : /3=0.5, M*=1OM . 



b : /3=0.5, Af*=2OAf . c : /3=0.5, M*=3OM . d : (3=1, 
M*=1OM . e : /3=1, M*=2OM . f : 0=1, M*=3OM . The 
+ symbols correspond to the most probable orbital inclination, 
confidence contours are drawn at the 25, 68 and 90% confidence 
levels. The dotted line is the "eclipse limit" : points above this line 
are excluded as it would imply an eclipse of the neutron star by 
its companion, which is not observed. 



2 2.2 

r = a + s 



2ascos(ip) 



(5) 



where ip is the angle subtended at the neutron star be- 
tween the radial and the observer directions, and a is the 
orbital separation. Finally, for a circular orbit with orbital 
inclination i, the angle ip is related to the inclination angle 
and the orbital phase <j> by 



typically lies in the range 1000-1500 km s 1 for a Bl star 
|Gathier et alj 1981). For ea rly-type stars, j3 is in the range 
0.7-1.2 Ce.g. lPuls et al.|[l996l ). We have Vo < Voo, so for our 
purposes we can take: 



cos(tp) = —sini cos(cj>(t) — 0(rgo)) 



(6) 



where rgo is the instant of the superior conjunction. 
Since no eclipse is visible in the system's light curve, we can 
restrict the inclination angle to 



V(r) ~ V a 



(2) 



i < arccosl — ) 
a 



(7) 



10 



L. Prat, J. Rodriguez, D.C. Hannikainen and S.E. Shaw 

Strong Absorption 



Low density 
highly ionized 
gas tail 



Phase 0. 75 



Neutron Star 



Soft Excess Visibility 




Phase 0.25 



Phase 




Weak Absorption 

Figure 7. Diagram of IGR J19140+0951 as it could be seen from the Earth. The orbital inclination of the 
system is taken to be ~ 65° . 



Table 4. Typical parameter values for a B1I star. 



Parameter 




value 


R* 


10 


- 30 R Q 


M< 


10 


- 30 Mq 


Ah 


10- 




Voc 


ss 1000 km s" 1 b 


p 




0.5-1 



a lLamersI <198ll) 

b lHaberl et all dl989t) 

4-3.2 Parameter fitting 

Our model needs four parameters, R+, M*, f3 and the ratio 
M*/V x , and computes the most probable orbital inclination 
i. We first verified the consistency of the ephemeris. We fixed 
every parameter to the typical values listed in Table [4] and 
let the program adjust the maximum of the absorption for 
several i. We obtained maxima between phases 0.49 and 0.5, 
thus confirming the value of rgo. Then, fixing the maximum 
to phase 0.5, we left free the inclination i and the stellar 
radius R t , set M* to 10, 20 and 30 Mq, f3 to 0.5 and 1, and 
the ratio M+/V00 normalized to the data. 

The results are shown in Fig. [6] The outer contour lines 
correspond to the 90% confidence levels. The two lines at 
2 Rq and 45 Rq are artifacts caused by the artificial con- 
straint of R+ between these two values. For the most prob- 
able stellar radius of 21 Rq, the lower inclination limit is 



constrained between 38° and 43°, with \ 2 being at a min- 
imum in the range 60-67°. Fig. left, shows the best-fit 
model against the experimental normalised data, with a 
good agreement. 

We can also determine the ratio K = -; — * „ . Using 
the values given in Table |4j the predicted value of K is Kq 
= 5.0 1 34 atoms cm . An estimate of the experimental 
value of K can be obtained using the best model, whose 
parameters are M* = 20M Q , R* = 21R Q and (3 = 0.5. 
We let the program adjust the normalisation of the model, 
and obtain K — 0.045 Kq (Fig. right). The observations 
predict Voo in the r ange 1000-1500 km s _1 for a B1I star 
|Gathier et al.lll98lh . which give M* in the range 4-7x 10~ 8 
MQ/year. This is consistent with the expected value, and 
retrospectively corroborates the model. However, the wind 
model used in this analysis is really simple, it does not take 
into account e.g. a clumping of the wind, since our data 
are not sufficiently precise to allow a refinement the model. 
Therefore our numerical results should remain mostly in- 
dicative of the evolution of the system. 



5 CONCLUSIONS: SIMPLE 

REPRESENTATION OF THE SYSTEM 

Fig.[7]summarizes schematically our main results. We found 
that photoelectric absorption highly correlates with the or- 
bital phase of the system. Using a simple stellar wind model 
we found a rather high orbital inclination, ~ 65°. We have 
detected a soft excess in some observations, just before the 
superior conjunction, in the area indicated on the figure. 
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This may be explained by a cloud of highly ionised gas sur- 
rounding the neutron star. Because of its ionisation, this 
gas is less accelerated by the stellar wind and, while being 
gradually overtaken by the wind, tends to form a "tail" . For 
typical parameters for the secondary star, the neutron star 
velocity Vns on its orbit is ~100 km s _1 and the wind ve- 
locity at the position of the orbit is ~0.8 Voo- Therefore, 
as the diagram shows, the angle between the "tail" and a 
vector normal to the orbit is 8 = arcsin( ( ^£!f ) « 10°. This 
tail could scatter the hard X-ray emission from the compact 
object and thus explain the soft excess feature. However, 
these results remain qualitative, as a precise determination 
of the geometry of the stellar wind would require a better 
sensitivity in the soft X-ray range, which could be achieved 
for instance with the use of the XMM or Chandra satellites. 

The study of X-ray binaries is challenging since it is 
often difficult if not impossible to identify their visible and 
infrared counterparts. Even if an infrared counterpart were 
observed, the distance to some systems prohibits the mea- 
surement of the orbital characteristics. Our study shows that 
X-ray observations can overcome these limitations and pro- 
duce very precise inferences. The RXTE and INTEGRAL 
observations of IGR J19140+0951 have led to good measure- 
ments of the orbital period of the system and constraints on 
its inclination angle. 

Moreover, we can use the compact object to probe 
the stellar wind of the companion. In the case of 
IGR J19140+0951, we diagnosed the type of the compan- 
ion (supergiant O or B), the wind density and its structure 
around the neutron star. More precise observations could 
lead to constraints on the mass and radius of the companion, 
and better constraints on the stellar wind. This allows the 
study of new INTEGRAL sources, either distant or highly 
absorbed, and ultimately the determination of new useful 
data for X-ray binary evolution scenarios. 
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